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Abstract 

The  linear  stability  of  circular  Couette  flow  between  concentric  infinite  cylinders 
is  considered  for  the  case  that  the  stationary  outer  cylinder  is  a crystal-melt  interface 
rather  than  a rigid  surface.  A radial  temperature  difference  is  maintained  across  the 
liquid  gap,  and  equations  for  heat  transport  in  the  crystal  and  melt  phases  are  included 
to  extend  the  ordinary  formulation  of  this  problem.  The  stability  of  this  two-phase 
system  depends  on  the  Prandtl  number.  For  small  Prandtl  number  the  linear  stability 
of  the  two-phase  system  is  given  by  the  classical  results  for  a rigid-walled  system.  For 
increasing  values  of  the  Prandtl  number,  convective  heat  transport  becomes  significant 
and  the  system  becomes  increasingly  less  stable.  Previous  results  in  a narrow-gap 
approximation  are  extended  to  the  case  of  a finite  gap,  and  both  axisymmetric  and 
non-axisymmetric  disturbance  modes  are  considered.  The  two-phase  system  becomes 
less  stable  as  the  finite  gap  tends  to  the  narrow-gap  limit.  The  two-phase  system 
is  more  stable  to  non-axisymmetric  modes  with  azimuthal  wavenumber  n — 1;  the 
stability  of  these  n — 1 modes  is  sensitive  to  the  latent  heat  of  fusion. 
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1.  Introduction 


For  crystal  growth  from  the  melt,  the  quality  of  crystal  produced  depends  strongly  on  the 
shape  of  the  crystal-melt  interface  during  the  growth  and  on  the  transport  of  heat  and 
mass  to  the  interface  from  the  nearby  fluid.  In  particular,  instabilities  of  the  interface  [1] 
and  hydrodynamic  instabilities  [2]  can  promote  undesirable  inhomogeneities  in  the  growing 
crystal.  Considerable  effort  goes  into  the  modelling  of  crystal  growth  techniques  in  order  to 
understand  and  control  the  process.  Since  the  theoretical  problem  involves  heat  and  mass 
transfer  via  nonlinear  convection-diffusion  equations  in  complex  geometries,  coupled  with  the 
free  boundary  representing  the  crystal-melt  interface,  many  studies  have  considered  simple 
geometries  for  which  there  are  closed-form  solutions  of  the  Navier-Stokes  equations. 

A crystal-melt  interface  bounding  a thermally-supercooled  [3]  or  “constitutionally-supercooled” 
[4]  fluid  is  subject  to  morphological  instability  [5],  and  one  line  of  inquiry  is  to  examine  the 
effect  of  various  types  of  flow  on  this  interfacial  instability  [6].  For  example,  the  effect  of 
plane  Couette  flow  [7,  8],  rotating  disk  flow  [9],  plane  stagnation  flow  [10,  11],  thermosolu- 
tal  convection  [12,  13,  14],  and  the  asymptotic  suction  profile  [15],  on  the  instability  of  a 
constitutionally-supercooled  interface  have  been  investigated.  Another  possibility  is  to  ex- 
amine the  characteristics  of  classical  hydrodynamic  instabilities  when  a bounding  surface  is 
a crystal-melt  interface  rather  than  a rigid  wall.  Such  instabilities  as  Benard  convection  [16], 
double-diffusive  convection  [12,  17,  18],  Poiseuille  flow  [19],  the  asymptotic  suction  profile 
[15],  and  thermal  and  shear  instabilities  in  buoyant  flow  due  to  lateral  heating  [19,  20,  21] 
have  all  been  examined  in  this  light.  In  general,  the  problem  of  predicting  the  extent  of 
the  coupled  interaction  between  the  crystal-melt  interface  and  the  flow  field  is  not  well- 
understood. 

In  1923  G.I.  Taylor  considered  the  stability  of  the  steady  Couette  flow  between  rotating 
cylinders  and  discovered  a transition  to  steady  axisymmetric  toroidal  rolls  (“Taylor  vor- 
tices”), obtaining  agreement  between  the  measured  values  of  the  critical  rotation  rates  and 
his  theoretical  predictions  [23,  24].  If  there  is  a radial  temperature  difference  across  the 
cylindrical  gap,  the  transport  of  heat  is  altered  by  the  formation  of  Taylor  vortices,  since 
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convective  heat  transport  then  augments  the  radial  diffusive  transport.  If  one  of  the  bound- 
aries is  a crystal-melt  interface,  this  change  in  heat  transport  can  modify  the  interface  shape. 
More  importantly,  this  change  in  morphology  can  alter  the  conditions  for  the  onset  of  the 
Taylor- vortex  flow. 

In  a previous  paper  we  considered  the  effect  of  a crystal-melt  interface  on  the  Taylor- 
vortex  flow  in  the  narrow-gap  approximation  [22].  In  this  note  we  extend  the  investigation 
to  include  the  effects  of  a finite  gap  width  and  non-axisymmetric  disturbances. 

2.  Theory 

The  geometry  is  shown  in  Fig.  1;  we  consider  the  linear  stability  of  an  idealized  system  with 
infinite  aspect  ratio  (— oo  < 2 < oo  in  cylindrical  variables).  The  crystal- melt  interface  is 
located  at  r = h(z,(f>,t ),  where  in  the  unperturbed  base  state  h{z,(j),t ) = R 2 is  constant. 
The  heated  inner  cylinder  of  radius  R\  rotates  at  angular  velocity  fix,  setting  up  a base  state 
consisting  of  a circular  Couette  flow  in  the  annular  melt  occupying  the  region  R\  < r < R2. 
The  melt  is  bounded  on  the  outside  by  a stationary  annulus  of  crystal  occupying  the  region 
R2  < r < R2  -f  Ls-  Both  the  heated  inner  cylinder  at  r — R\  and  the  cooling  bath  at 
r = R2  + Ls  are  each  assumed  to  be  isothermal,  as  is  the  crystal- melt  interface,  which  is  at 
the  equilibrium  melting  temperature.  The  densities  of  the  melt  and  crystal  are  assumed  to 
be  the  same,  and  the  radial  temperature  difference  across  the  melt  is  assumed  to  be  small 
enough  that  any  buoyancy  effects  [25]  can  be  neglected. 

For  details  we  refer  to  Ref.  [22].  Briefly,  the  nonlinear  governing  equations  are  the 
incompressible  Navier-Stokes  equations  for  the  velocity,  it,  and  the  pressure,  p,  in  the  melt 
region  Rx  < r < h(z,<f>,t), 

V • it  = 0, 

du  . 1 ~ 

— — b (it  • V)it  H — Vp  = t^V2it, 

ot  p 

and  the  convection-diffusion  equation  for  the  temperature  in  the  melt,  T, 


(la) 

(lb) 


and  the  diffusion  equation  for  the  temperature  in  the  crystal,  Ts, 


dTs 

dt 


KsV2Ts. 


(3) 


Here  p is  the  density,  v is  the  kinematic  viscosity,  and  k and  k s are  the  thermal  diffusivities 
of  the  melt  and  crystal,  respectively.  At  the  crystal-melt  interface  r = h(z,<f),t)  and  the 
inner  cylinder  r = R\  the  flow  satisfies  no  slip  and  no  normal  flow,  u = 0.  The  temperatures 
are  specified  at  the  inner  cylinder  and  the  outer  edge  of  the  crystal.  The  thermal  boundary 
conditions  at  the  interface  are 


,.dT  . dTs  _ r 

™ r\  r\  LyVn 

on  on 

(4a) 

T = Ts  — Tm  + TmTK , 

(4b) 

where  d/dn  is  the  normal  derivative,  k and  ks  are  the  thermal  conductivities  of  the  liquid 
and  crystal,  respectively,  Tm  is  the  bulk  melting  point  of  the  crystal,  Ly  is  the  latent  heat 
of  fusion  per  unit  volume  of  crystal,  vn  is  the  normal  velocity  of  the  interface,  and  T is  a 
capillary  length.  Eq.  (4a)  expresses  the  conservation  of  energy:  the  latent  heat  released  by 
the  phase  transition  is  conducted  into  the  crystal  and  melt.  Eq.  (4b)  is  the  Gibbs-Thomson 
equation  which  relates  the  interface  temperature  to  the  mean  curvature  of  the  interface,  1C. 

We  consider  the  stability  of  the  one-dimensional  base  state  given  by  h(z,(f>,t)  = R 2, 
u = w = 0,  v — W°)(r),  p = pf°f(r),  T — T^°^(r),  and  Ts  = Tg°\r),  where  u,  v , and  w are 
the  radial,  azimuthal,  and  axial  velocity  components,  respectively, 


y(°)(r)  = rili  — 

U (Rl/Rl- 1) 

r<°>(r)  = Tm  + rm-£-  + Ar'0^1 

R2  \og{R\l  R2) 


and 


r<»  = rm  + rmi-+(yAr; 


log  (r/R2 


log  (RJRa)’ 

the  base  pressure  gradient  balances  the  centripetal  acceleration  term  p(V^)2/r. 


(5) 

(6) 

(7) 
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To  determine  the  linear  stability  of  the  base  state,  we  assume  the  dependence  of  the 
perturbed  quantities  takes  the  form  f(r)  exp (crt  + in<t>  + m),  and  write 


“ 1 

f ° ) 

^ u(r)  ^ 

V 

y(°)(r) 

v(r) 

w 

0 

w(r) 

p 

= 

p(°)(r) 

+ 

P(r) 

exp(o-f  + in<f)  + iaz), 

T 

r<°)(r) 

T(r) 

Ts 

rfV) 

Ts(r) 

h ) 

V R‘  / 

{ * ) 

where  the  small  quantities  w,  v,  w,  p , T,  Ts  and  6 are  complex.  The  exponential  growth 
rate  a is  complex  as  well,  and  the  base  state  is  linearly  stable  if  ay,  the  real  part  of  <j,  is 
negative  for  all  wavenumbers  a and  n.  A marginal  disturbance  (ay  = 0)  is  stationary  if  o’,-, 
the  imaginary  part  of  a,  vanishes. 

The  equations  and  boundary  conditions  are  linearized  about  the  base  state  and  the 
resulting  eigenvalue  problem  describing  the  hnear  stability  of  the  system  is  solved  numerically 
using  the  same  approach  described  in  [22].  The  dimensionless  parameters  which  enter  the 
problem  are  the  Taylor  number 


the  Prandtl  number  Pr  = z///c,  the  ratio  Ps  = v/ks,  the  dimensionless  latent  heat  of  fusion 
jC  = vLv /(fcAT),  the  dimensionless  surface  energy  7 = TmT  / ([R2  — i?i]AT),  the  ratio  of 
thermal  conductivities  q — ks/k,  the  radius  ratio  7 = Ri/R2,  and  the  dimensionless  crystal 
width  Lc  = Ls/[R2  ~ #1]. 


3.  Numerical  Results 

Our  principal  concern  is  the  interpretation  of  recent  experiments  with  succinonitrile  (SCN) 
with  Pr  = 22.8,  for  the  case  of  a finite  gap  7 < 1.  For  centimeter-sized  systems  with 
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moderate  temperature  differences  across  the  melt  gap,  the  parameter  7 is  very  small  and 
can  be  neglected.  The  thermal  properties  of  the  crystal  and  melt  are  very  similar  for  SCN, 
and  we  will  take  q — 1 and  Ps  = 22.8.  We  take  Lc  = 1;  the  effect  of  varying  Lc  is  considered 
in  Ref.  [22]  for  the  narrow  gap  limit.  We  will  present  our  numerical  results  in  dimensionless 
form,  using  the  length  scale  [R2  — i?i],  the  time  scale  [R2  — Ri]2/v,  the  temperature  scale 
AT,  and  the  fluid  velocity  scale  iT  1 [T?2  — i?i].  However,  we  will  retain  the  previous  notation 
for  all  quantities,  which  will  henceforth  be  dimensionless. 

3.1.  Axisymmetric  modes 

For  axisymmetric  disturbances  to  the  two-phase  system  we  find  that,  as  in  the  rigid-walled 
system,  the  onset  of  instability  is  stationary,  with  a,  = 0.  In  this  case  the  marginal  Taylor 
numbers  do  not  depend  on  C,  Ps,  or  q.  In  Fig.  2 we  show  marginal  stability  curves  for  a 
fixed  radius  ratio  q — 0.5  for  various  values  of  the  Prandtl  number.  For  PT  = 0 the  thermal 
transport  is  diffusion-controlled  and  convection  has  no  effect  on  the  thermal  field.  The 
isotherms,  and  hence  the  crystal- melt  interface,  are  therefore  cylindrical,  as  dictated  by  the 
isothermal  boundary  conditions  at  r = Ri  and  r = R2  + L.  The  results  for  Pr  = 0 therefore 
are  identical  to  the  results  for  a rigid- walled  system,  with  a critical  Taylor  number  Ta  = 6,  200 
corresponding  to  a dimensionless  wavenumber  a = 3.2.  As  Pr  increases,  convective  transport 
of  heat  becomes  significant,  and  the  interface  responds  to  the  influence  of  the  flow  field,  which 
allows  the  instability  to  set  in  sooner:  for  Pr  = 22.8,  the  critical  disturbance  has  Ta  = 1,  060 
and  a — 1.7. 

In  Figs.  3 and  4 we  show  the  linear  eigenmodes  for  Pr  = 0.01  and  Pr  = 22.8  in  detail, 
together  with  contour  plots  of  the  flow  field  and  temperature  in  the  melt.  On  the  left  of 
each  plot  we  show  the  dimensionless  angular  velocity  of  the  base  state  H(r)  = W0)(r)/r, 
and  the  dimensionless  eigenfunctions  u(r),  e(r),  and  T(r),  normalized  so  that  the  maximum 
value  of  u is  one.  The  contour  plots  show  a cross-section  of  the  melt  and  crystal,  with  the 
cross-hatched  region  corresponding  to  the  crystal;  the  plots  are  obtained  by  adding  to  the 
base  state  a small  multiple  (0.003)  of  the  normalized  eigenmode,  in  order  to  highlight  the 
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interface  deformation.  For  Pr  — 0.01  with  a — 3.162  the  effect  of  the  crystal-melt  interface 
on  the  classical  Taylor-vortex  flow  is  small,  and  the  flow  consists  of  counter-rotating  cells 
of  roughly  unit  aspect  ratio.  The  isotherms  are  undisturbed  by  the  flow,  and  show  the 
logarithmic  spacing  of  the  base  state.  For  Pr  = 22.8  with  a = 1.697  (Fig.  4)  the  perturbed 
azimuthal  velocity  and  temperature  fields  are  larger,  and  the  interface  and  isotherms  are 
distorted  by  the  flow  field.  The  counter-rotating  cells  have  a rectangular  cross-section.  The 
interface  melts  back  where  the  perturbed  flow  transports  additional  heat  to  the  interface, 
and  protrudes  into  the  liquid  where  the  flow  carries  heat  away  from  the  interface.  Note  also 
that  v is  much  larger  near  the  interface  for  Pr  = 22.8  than  it  is  in  the  rigid  case  for  Pr  = 0. 

In  Fig.  5 we  show  marginal  stability  curves  for  Pr  = 22.8  for  various  values  of  the  radius 
ratio  77.  As  is  the  case  for  the  rigid- walled  system,  the  effect  of  decreasing  7/  from  the  narrow- 
gap  limit  77  = 1 is  to  stabilize  the  system.  As  can  be  seen  from  Table  I,  the  amount  of 
destabilization  relative  to  the  rigid- walled  system  remains  about  the  same  over  this  range  of 
77,  amounting  to  a five-fold  to  seven-fold  decrease  in  the  critical  Taylor  number;  the  critical 
wavenumber  is  not  sensitive  to  77  over  this  range. 

3.2.  Non-axisymmetric  modes 

We  next  consider  the  behavior  of  non-axisymmetric  disturbances  with  n ^ 0.  For  a rigid- 
walled  system  in  the  narrow  gap  limit,  the  most  dangerous  mode  is  axisymmetric,  but  there 
are  also  nearby  non-axisymmetric  modes  at  only  slightly  higher  Taylor  numbers  [26].  In 
contrast  to  the  axisymmetric  case,  which  has  a stationary  onset  of  instability  (a,  = 0),  the 
non-axisymmetric  modes  have  non-zero  values  of  <7,-,  and  the  disturbance  propagates  as  a 
travelling  wave.  Thus,  whereas  the  critical  Taylor  number  for  the  axisymmetric  mode  is 
insensitive  to  the  latent  heat  in  the  two-phase  problem,  for  the  non-axisymmetric  mode 
energy  must  be  continually  supplied  locally  to  transform  crystal  to  melt  to  sustain  the 
progressing  wave.  Consequently  the  critical  Taylor  numbers  for  the  non-axisymmetric  modes 
are  found  to  depend  on  the  value  of  the  latent  heat. 

The  dimensionless  parameter  C = uLy / (kAT)  is  a measure  of  the  amount  of  energy 
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required  to  transform  crystal  into  melt.  Either  a small  value  of  the  latent  heat  Ly  or  a large 
temperature  difference  AT  across  the  melt  give  rise  to  a small  value  of  £,  meaning  that  it  is 
relatively  easy  to  melt  the  crystal.  If  the  latent  heat  is  large,  or  the  temperature  difference 
across  the  melt  is  small,  it  is  difficult  to  provide  enough  heat  to  melt  the  crystal.  Evaluating 
the  material  constants  for  SCN,  the  dimensionless  latent  heat  can  be  written  in  the  form 
C = 557. 3/ AT,  for  AT  given  in  Kelvins.  As  shown  in  Fig.  6 for  77  = 0.99,  for  n — 1 modes 
the  two-phase  system  is  least  stable  for  C = 0,  becoming  more  stable  with  increasing  C.  The 
marginal  Taylor  number  for  the  n = 1 modes  for  C = 0 is  Ta  = 1150,  which  is  significantly 
higher  than  the  n = 0 value,  Ta  — 670.  As  Ta  is  increased,  the  neutral  stability  curves  begins 
to  pinch  off,  as  shown  for  C = 50,  and  eventually  forms  two  branches,  as  shown  for  C = 350. 
The  upper  of  the  two  branches  tends  to  the  marginal  curve  corresponding  to  the  rigid- walled 
system  for  large  values  of  £,  and  the  interface  deflection  for  the  normalized  eigenmode  tends 
to  zero  as  £-1 . The  interface  deflection  for  modes  on  the  lower  closed  loop  of  marginal  states 
tends  to  a non-zero  limit  for  large  £,  and  the  product  Ccrt  tends  to  a constant.  We  have 
not  investigated  further  the  nature  of  this  lower  branch,  as  in  any  event  the  most  dangerous 
mode  for  the  two-phase  system  corresponds  to  an  axisymmetric  mode  for  these  parameter 
values. 

For  non-axisymmetric  modes,  then,  the  destabilizing  effect  of  convective  heat  transport 
to  the  interface  at  large  Prandtl  number  can  be  somewhat  offset  by  large  values  of  C.  In 
Table  II  we  compare  the  marginal  Taylor  numbers  for  n = 0 and  n = 1 modes  in  the  rigid- 
walled  system  ( Pr  = 0)  and  the  two-phase  system  (Tr  = 22.8)  for  radius  ratios  77  near  unity. 
The  comparison  is  made  with  C = 0,  which  is  the  least  stable  case  for  the  n — 1 mode.  For 
Pr  = 22.8,  the  n = 1 mode  is  substantially  more  stable  than  the  n = 0 mode,  unlike  the 
rigid-walled  case.  The  n = 1 mode  is  also  much  more  stabilized  for  decreasing  77  than  are 
the  other  modes  in  the  systems. 


-8- 


4.  Discussion 


The  destabilization  of  the  two-phase  system  relative  to  the  rigid-walled  system  for  large 
Prandtl  numbers  that  was  found  for  the  case  of  a narrow  gap  [22]  also  occurs  for  a finite 
gap  system,  although  the  relative  destabilization  is  largest  for  the  narrow  gap  case.  As  for 
the  case  of  a rigid- walled  system,  for  steady  rotation  of  the  inner  cylinder  with  a stationary 
outer  cylinder,  the  onset  of  instability  is  due  to  a stationary,  axisymmetric  mode.  The  non- 
axisymmetric  modes  with  azimuthal  wavenumber  n — 1 that  are  nearby  the  onset  for  the 
rigid-walled  system  are  found  to  be  relatively  more  stable  for  the  two-phase  system. 

A possible  extension  of  this  work  would  include  the  effects  of  buoyancy-driven  flow  caused 
by  the  radial  temperature  gradient  in  a vertical  alignment  of  the  system.  As  in  the  rigid- 
walled  case  [27],  a one-dimensional  solution  may  be  found  for  the  limiting  case  of  an  infinite 
aspect  ratio,  and  a tractable  linear  stability  problem  may  be  formulated. 

It  is  also  of  interest  to  consider  the  case  of  a time- varying  torsional  oscillation  of  the  sys- 
tem with  zero  mean  oscillation,  and  preliminary  experimental  and  theoretical  work  on  this 
problem  is  in  progress;  the  theoretical  treatment  is  based  on  Floquet  theory.  Preliminary 
experimental  results  for  the  two-phase  system  using  succinonitrile  [22]  suggest  an  axisym- 
metric mode  of  instability,  whereas  the  numerical  results  of  Carmi  and  Tustaniwskyj  [28]  for 
an  oscillating  rigid- walled  system  indicate  a non-axisymmetric  (n  — 1)  mode  of  instability, 
with  nearby  axisymmtric  modes  that  are  only  slightly  more  stable. 
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Table  I.  Comparison  of  Rigid  and  Two-Phase  Systems 


Pr 

rj  — 0.99 

Ta  a 

r i = 0.50 

Ta  a 

tj  = 0.30 

Ta  a 

t)  = 0.20 

Ta  a 

0.0 

3,412.63  3.127 

6,199.16  3.162 

11,548.16  3.215 

20,715.77  3.263 

22.8 

670.14  1.749 

1,059.17  1.697 

1,834.04  1.671 

3,144.69  1.652 

Table  II.  Comparison  of  Axisymmetric  and  Non-Axisymmetric  Modes 


Pr 

T]  = 

n — 0 

0.99 

n = 1 

V = 

n = 0 

0.97 

n — 1 

V = 

n = 0 

0.95 

n = 1 

T 

3,412.63 

3,416.93 

3,459.89 

3,473.51 

3,509.71 

3,533.70 

0.0 

a 

3.127 

3.128 

3.127 

3.131 

3.127 

3.135 

cr. 

0.0 

-2.180 

0.0 

-3.828 

0.0 

-5.009 

T 

670.14 

1,150.46 

676.07 

1,692.84 

682.39 

2,081.06 

22.8 

a 

1.749 

2.322 

1.747 

2.636 

1.746 

2.787 

<7 1 

0.0 

-0.398 

0.0 

-0.819 

0.0 

-1.236 
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Figure  Captions 


Figure  1.  Schematic  of  the  case  of  a stationary  crystalline  annulus  (c)  surrounding  an 
annular  melt  (^)  in  contact  with  an  inner  cylinder  rotating  at  angular  velocity 

Figure  2.  Marginal  Taylor  numbers  versus  the  axial  wavenumber  a for  axisymmetric  modes 
with  Prandtl  numbers  Pr  = 0.0,  5.0,  22.8,  and  100.0. 

Figure  3.  The  linear  eigenmode  for  Pr  = 0.01  and  77  = 0.5.  0(r),  u(r),  5(r),  and  T(r ) are 
the  dimensionless  angular  velocity  of  the  base  state,  the  radial  part  of  the  perturbed  dimen- 
sionless radial  velocity,  the  radial  part  of  the  perturbed  dimensionless  azimuthal  velocity, 
and  the  radial  part  of  the  perturbed  dimensionless  temperature  in  the  melt,  respectively. 
Streamlines  of  the  flow  are  shown  in  the  contour  plot  in  the  middle,  with  the  cross-hatched 
region  corresponding  to  the  crystal.  Isotherms  are  shown  in  the  contour  plot  on  the  right. 

Figure  4.  The  linear  eigenmode  for  PT  = 22.8  and  77  = 0.5.  As  in  Figure  3,  in  the  contour 
plots  the  linear  perturbation,  normalized  so  that  max{u(r)}  = 1,  is  added  to  the  base  state 
with  an  amplitude  of  0.003  to  show  the  interface  deformation. 

Figure  5.  Marginal  Taylor  numbers  for  axisymmetric  modes  with  Pr  = 22.8,  versus  the 
axial  wavenumber  a for  radius  ratios  77  = R1/R2  of  0.990,  0.50,  0.30,  and  0.20. 

Figure  6.  Marginal  Taylor  numbers  for  non-axisymmetric  modes,  versus  the  axial  wavenum- 
ber a for  values  of  the  dimensionless  latent  heat  C = 0,  50,  and  350.  There  are  two  branches 
of  the  marginal  stability  curve  for  C = 350;  the  upper  branch  is  nearly  identical  to  that  for 
the  rigid- walled  system. 
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